Low-Energy Structures in Strong Field Ionization Revealed by Quantum Orbits 
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Experiments on atoms in intense laser pulses and the corresponding exact ab initio solutions of the time- 
dependent Schrodinger equation (TDSE) yield photoelectron spectra with low-energy features that are not re- 
produced by the otherwise successful work horse of strong field laser physics: the "strong field approximation" 
(SFA). In the semi-classical limit, the SFA possesses an appealing interpretation in terms of interfering quantum 
trajectories. It is shown that a conceptually simple extension towards the inclusion of Coulomb effects yields 
very good agreement with exact TDSE results. Moreover, the Coulomb quantum orbits allow for a physically 
intuitive interpretation and detailed analysis of all low-energy features in the semi-classical regime, in particular 
the recently discovered "low-energy structure" [C.I. Blaga et al, Nature Physics 5, 335 (2009) and W. Quan et 
al, Phys. Rev. Lett. 103, 093001 (2009)]. 



PACS numbers: 32.80.Rm, 32.80.Wr, 34.80.Qb 

The development of analytical and numerical methods ca- 
pable of treating strongly-driven quantum systems is of great 
interest in many areas of physics. By "strongly-driven" we un- 
derstand that conventional time-dependent perturbation the- 
ory is not applicable. A prime example for such a system is 
an atom in an intense laser field. The force on valence elec- 
trons due to the electric field of the electromagnetic wave de- 
livered by present-day intense lasers can easily compete with 
the binding force. As a consequence, the photoelectron spec- 
tra may show strong nonperturbative features such as plateaus 
and cut-offs [ 1 ], instead of a simple exponential decrease with 
the number of absorbed photons, as expected from perturba- 
tion theory. Recently, an "ionization surprise" (2|] at wave- 
length A = 2 p,m and intensity I = 80-150 TW/cm 2 , the so- 
called "low-energy structure" (LES) JHQ], has been reported. 
The LES is a strong but narrow enhancement of the differen- 
tial ionization probability along the polarization direction of 
the laser at low energies. This result was so astonishing not 
only because it is unpredicted by the "strong field approxi- 
mation" (SFA) but also because it is observed in a regime 
where matters were actually expected to simplify. In fact, if 
the number of photons N of energy fiu) required to overcome 
the ionization potential I p is large, N = I p /huj 3> 1, and 
the time the electron needs to tunnel through the Coulomb- 
barrier is small compared to a laser period, i.e., the Keldysh 
parameter 7 = */ I p /2U P with U p the ponderomotive poten- 
tial, is small, a quasi-static tunneling theory appears to apply 
|0]. As the tunneling ionization rate in a static electric field is 
a smooth, featureless function of the final momenta pit and p± 
parallel and perpendicular to the electric field, respectively, 
no LES has been expected. In the present Letter we reveal the 
origin of the LES using our trajectory-based Coulomb-SFA 
(TC-SFA). The fact that the TC-SFA allows recourse to tra- 
jectories provides an unprecedented insight into the origin of 
any spectral feature of interest, as constructive or destructive 



interference of trajectories or the Coulomb-focusing of them 
10] can be analyzed in all details. 

The SFA is a widely and successfully used theoretical ap- 
proach to tackle strong field photo detachment. In its sim- 
plest form it accounts only for the so-called "direct" elec- 
trons, which are bound to the atom up to the detachment 
time to and at later times just move in the laser field with- 
out any interaction with the parent atom anymore. Quantum 
mechanically, such a free electron with drift-momentum p in 
a laser field defined by the vector potential A(t) is described 
in length gauge by a Gordon- Volkov state JH] |\l/p GV ' ) (t)) — 
e -iS p (t)| p + A ^ Sp ^ = + A(t')] 2 /2dt' (atomic 
units are used unless noted otherwise). The SFA-transition 
matrix element for a final asymptotic momentum p (at the de- 
tector) for a laser pulse E(t) — —dtA(i) in dipole approxima- 
tion, polarized along the z-axis, and lasting from t = until 
t = T p then reads 



M p SFA > = -i 



{¥ p GV \t)\zE(t)\y a (t))dt (1) 



where 1*0(0) = e i/pt |*o>, and |<J/ ) is the initial bound state 
of the electron. Alternatively, we may use the matrix element 
^yKSFA) _ g-iSpfTpj-UpTp j^£SFA)^ as ^ gj ves j ne same pho- 
toelectron spectrum, dw(p)/d 3 p = |M p SFA ^| 2 = |M p SFA ^| 2 . 
The neglect of any further interaction between electron and 
binding potential once the electron is ejected is well-justified 
for short-range potentials as, e.g., in the photodetachment 
from negative ions |^[l0[]. However, the agreement with pho- 
toelectron spectra calculated ab initio by solving the TDSE is 
in general poor for the case of long-range binding potentials, 
such as in the more common ionization of neutral atoms or 
positive ions 111 [Till . For illustration, we show in Fig. [T] SFA 
and TDSE photoelectron momentum spectra for a linearly po- 
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FIG. 1: Logarithmically scaled photoelectron momentum distribu- 
tions in the p z pa;-plane (py = p z ) calculated using (a) plain SFA, (b) 
TC-SFA, and (c) numerical solution of the TDSE. Panel (d) shows 
the linearly scaled, angle-resolved energy spectrum in backward- 
direction with the LES around 1 eV clearly visible. The laser and 
atomic parameters are given in the text. Each spectrum was normal- 
ized to its maximum value. The small empty semi-circle around the 
origin in the TDSE result (c) is due to the mapping of the angle- 
resolved energy spectra (as calculated by Qprop 1 12]) to the Carte- 
sian p z p x -momentum plane. 



larized n — 3-cycle pulse of the form 
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for t £ [0, 2mr/cj], and zero otherwise. The peak field 
strength E = 0.0534 corresponds to 100TW/cm 2 , and the 
laser frequency to = 0.0228 to A = 2 /jm. The binding 
potential used in the TDSE solver Qprop |[l2Tl was V(r) — 
-l/r - 17.0e~ 17A3r /r, which yields I p = 0.579 for the 
Is state and a grid spacing Ar — 0.2. The same I p and a 
Coulombic Is initial state was used in the SFA-calculation. 
Theparticular choice of I p was motivated by the experiment 
in ||3|,[l3|] where argon was used. Figure[TJshows several strik- 
ing discrepancies between SFA and TDSE: (i) the SFA pre- 
dicts a symmetric momentum distribution for a pulse of the 
form (O with an integer number of cycles n while the TDSE 
spectrum is strongly asymmetric 11 1 311 - (ii) The radial struc- 
tures lUll [lrjll around p = present in the TDSE spectrum 
are absent in the SFA result, (iii) Several side-lobes in the 
TDSE result are clearly visible for p z < but completely 
absent in the SFA. In Fig. [Tp we anticipate the spectrum ob- 
tained with our TC-SFA method, which is — as regards points 
(i)-(iii) above — in excellent agreement with the TDSE result. 
Only the probability along a clearly visible caustic is overes- 
timated. However, the classical caustic turns out to be related 
to the LES, as will be discussed below. 

The TDSE momentum spectrum in Fig. QJ does not al- 
low a clear identification of the LES. In order to compare the 
TDSE results with Refs. [Ht] angle-resolved energy spec- 
tra dw(£,9)/ sin9d9d£ = p\M p \ 2 (where £ = p 2 /2 and 
tan 9 — \p±]/p\\) were calculated. Figure [TJl clearly shows 
the LES for the directional energy spectrum anti-parallel (9 — 



180°) to the polarization axis. Because of the short pulse du- 
ration, no LES is visible for 9 = (not shown), which helps 
to identify unambiguously the origin of the LES in the follow- 
ing. Figure [TJl compares well with the experimental result in 

Fig. 3b of il. 

Let us now introduce the TC-SFA approach. The time in- 
tegral in (UJ can be evaluated semi-classically for sufficiently 
big N with the help of the saddle-point approximation leading 
to 



M^SFA) 



7 ,^pb' 



-iS. 



where S p i pS = S p i p (t ps ), S„i (t 



pips 



(3) 



A(t)} 2 /2 



I p }dt and C ps is a prefactor QjJ, L17D - Here, the sum is over all 
saddle-point times t ps fulfilling the stationary action equation 



dS. 



pip 



dt 



& -[p + A(t ps )f 



(4) 



For a given final momentum p there will, in general, be two 
complex saddle-point times t ps per laser cycle. The saddle- 
point action S p i pS can be separated into an integral from the 
complex saddle-point time t ps down to the real axis where 



5p/ p s 



to and then along the real axis from to to 
S^j s . In plain SFA both con- 
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tributions can be calculated analytically for a given (ana- 
lytic) vector potential. However, the TC-SFA requires the 
numerical evaluation of electron trajectories in the presence 
of both the binding potential and the laser field as well as 
the corresponding modified action along the real time axis 
Spf* s . Hence it is illustrative to first inspect such trajecto- 
ries for the plain SFA case. Consider the equations of mo- 
tion r = v = p + A(t) and p = 0. Clearly, p =const., 
and for a given saddle-point time t ps the real part of the ac- 
tion reads S^j s = J^ p {y 2 /2 + I p } dt, as required. If we 
actually want to solve for r(t) we need to fix the initial con- 
ditions for v and r. The initial condition for momentum (or 
velocity) follows from ©. For the initial position we choose 
Re r (tp S ) = 0, which implies r(trj) = a (to) — Rca:(tps) 
with a(t) = J A(t')dt'. In the tunneling-limit r(to) can 
be identified with the geometrical tunnel exit Zt(tn) given by 
-I p = E{t )z t {t ). 

In our TC-SFA the equation of motion p = is 
replaced by p = — r/|r| 3 , and simultaneously p is 
replaced by po in (|4j, leading to a modified saddle- 
point time t poS (and real part to = Ret poS ). The 
TC-SFA saddle-point action reads S„ T „ = 

r PO-lpS 

^ft s = &{[p + A{t)] 2 /2- 



I p } dt and 

S pJ p s = Si* - VK*) + h) d*- The action of the 

Coulomb-potential on the imaginary dynamics during the tun- 
neling process affects the ionization probability la, ll7[] but 
does hardly change the shape of the photoelectron spectra so 
that we neglect it in this work. However, note that due to the 
replacement of p by po and the Coulomb-modified saddle- 
point times nevertheless S^j} s ^ S pI s . 
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The TC-SFA "shooting"-method is implemented as fol- 
lows. Two loops run over poz £ [ — Pol ax ?Po" ax ] an ^ Pox £ 
[ — PQx &yL t Pox ax ]- Because of cylindrical symmetry about p|| = 
p z sampling of the p Qz p 0x -plane is sufficient. For each initial 
momentum Eq. (0]i is solved for t P()S using a complex-root- 
finding routine. Each t poS corresponds to a trajectory. Hence, 
there is another loop over all t PoS found. The complex part 
S s of the action for a given t poS is calculated analytically, 
neglecting Coulomb effects. The trajectory corresponding to 
a t PQS is calculated from to up to the end of the pulse ac- 
cording to the equations of motion z = p z + A, x = p x , 
p z = -z/(x 2 + z 2 ) 3 / 2 , p x = -x/(x 2 + z 2 ) 3 / 2 (all time- 
arguments suppressed) for the initial conditions x(to) = 0, 
z(t ) = a* (to) - Rea z (t poS ), p z {t Q ) = p 0z , p x (t Q ) = 
Po x using a Runge-Kutta solver. If the energy £(T P ) = 
p 2 (T p )/2 — \/r(T p ) is negative, the trajectory does not cor- 
respond to a free electron and thus does not contribute to the 
photoelectron spectrum If the energy is positive the 

asymptotic momentum p can be calculated from p(T p ) and 
r(T p ) using Kepler's laws [ 14], avoiding unnecessary explicit 
propagation up to large times. Finally, the result for the t pas 
under consideration is stored in a table of the form po z , po x , 
Pz, Px, tp s, Sp oI s , Once the loops over the initial mo- 
menta are completed, the table with the trajectory data can be 
post-processed. The trajectories are binned according to their 
asymptotic momentum, and the TC-SFA matrix element 

Mf c >=^C< P0S e- if W (5) 

s 

is calculated, where the sum is over all trajectories ending up 
in the final momentum bin centered at p. 

The TC-SFA result in Fig.QJ was obtained with ~ 2 x 10 8 
trajectories and shows the above introduced Coulomb-features 
(i)-(iii) of the TDSE calculation. In addition there is a caustic 
structure with a maximum at p c = (p c \\,Pc±) = (—0.22,0.1). 
We have checked that this maximum moves closer towards the 
polarization axis as the wavelength is increased (with all other 
laser and target parameters held constant). For shorter wave- 
lengths it fades away while becoming more ring-like. This is 
a strong indication that the LES of the full quantum TDSE 
simulations manifests itself as a caustic in the semi-classical 
TC-SFA approach. 

It was argued already in the accompanying article J2] 
that the LES is due to low-energy forward scattering at the 
Coulomb-potential. Certainly, only electrons emitted with 
high probability, i.e., when the absolute value of the electric 
field is high, can contribute to such a pronounced spectral fea- 
ture as the LES. The quantum trajectories of the TC-SFA re- 
sponsible for the caustic illustrate and confirm this view-point. 
For further analysis, we group all trajectories into four types. 
The first two types are the so-called "short" and "long tra- 
jectories" of plain SFA. The short trajectory (type I) fulfills 
ztPz > and p x po x > 0, meaning that the ejected electron 
moves from the tunnel exit directly towards the detector and 
has no close encounters with the ion. Coulomb effects on this 
type of trajectory are expected to be small. Type II is a "long 




FIG. 2: (color online), (a) Partial spectra due to trajectories I-IV. 
Panel (b) shows A(t) (solid) and E(t) (dotted) and illustrates why 
there is an interference structure for p z < in (a) for type-I trajec- 
tories but not for type II. (c) The dominant trajectories in the zx- 
plane up to the time T p which contribute to the final momentum p c , 
i.e., where the caustic structure is visible for the type-Ill spectrum in 
panel (a). 



trajectory" obeying z t p z < and p x po x > 0. Here, the elec- 
tron starts from the tunnel exit that points away from the de- 
tector but ends up with a parallel momentum in the opposite 
direction because of the drift it acquires from the laser field at 
the time of emission, i.e., in plain SFA p z — po z = —A(to). 
The lateral initial momentum po x of type-I and type-II tra- 
jectories is already in the same direction as the final momen- 
tum at the detector. This is always the case in plain SFA, as 
Px = Pox- Hence there are only type-I and type-II trajectories 
in plain SFA. While in Ref. [ 17] the type-I and type-II trajecto- 
ries existing already in plain SFA are Coulomb-corrected, the 
shooting-method introduced in the present Letter goes beyond 
a correction by taking into account qualitatively new types of 
trajectories generated by the binding potential. Type III is 
classified via z t p z < and p x po x < 0. The corresponding 
electrons start on the opposite tunnel exit, as type-II electrons 
do, but moreover have an initial lateral momentum pointing in 
the "wrong" direction. Only due to relatively close encoun- 
ters with the ion, the lateral momentum is ultimately reversed. 
Finally, type-IV trajectories obey z t p z > and p x po x < 0. 

In order to understand the origin of spectral features such 
as the LES it is useful to study the partial contributions of 
each class of trajectories to the total spectrum in Fig. [TJ). To 
that end we split the sum over all trajectories in (O into four 
terms, Mp C ' = Y^JLi Mup G \ where the lower index v in- 
dicates the trajectory type. Figure^ shows the partial spectra 
I -M^p ' \ 2 ■ The spectrum |M I < p rC ' 1 1 2 allows only trajectories of 
type I to interfere with each other. Figure |2j) illustrates why 
there is only an interference pattern for p z < but none for 
p z > 0: because of p z ~ — A(to) the vector potential A at 
the time of emission for p z < is positive. For a positive 
A there are two contributions (filled black circles in Fig. [2}}). 
The tunnel exit at those two times is at z < 0, so that in- 
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trix element including Coulomb effects. The method yields 
an unprecedented agreement with the full ab initio solutions 
of the time-dependent Schrodinger equation. The method al- 
lows to analyze any spectral feature in terms of interfering 
or Coulomb-focused quantum trajectories. In particular, the 
recently observed low-energy structure at long wavelengths 
was found to originate from low-energy forward scatter- 
ing, leading to caustic structures in semi-classical trajectory- 
calculations. 



FIG. 3: (color online). Emission times of type-Ill trajectories con- 
tributing to final momenta where the caustic structure is located. At 
p c the number of contributing trajectories is particularly high (steep 
slope). 

deed the two trajectories are both of type I. The open circle in 
Fig. indicates a type-II trajectory because z t > 0. There 
is just one dominant type-II trajectory, which explains that 
there is no interference pattern in the spectrum |Aifj I p C ' l | 2 for 
p z < in Fig. [2^. Figure [2J; shows the dominant trajectories 
in the zx-plane which contribute to the final momentum p c , 
i.e., where the caustic structure is visible for the type-Ill spec- 
trum in panel Fig. [2^. As is evident from Fig. [2}:, the type-Hi 
electrons are driven back close to the origin and slowed down 
by the laser field. 

We finally analyze how the caustic structure in the spectrum 

Mjjjr? 1 2 °f Fig. is formed. To that end we added all type- 
Ill contributions incoherently and obtained a spectrum very 
similar to \Mj^' | 2 of Fig. [2^. Hence, interference can be 
ruled out, and, indeed, the LES has been observed in classical 
ensemble calculations as well [4]. Our TC-SFA simulation 
reveals that the caustic structure is due to an enhancement of 
the number of trajectories of type III. This is shown explicitly 
in Fig. [3] where the emission times t of type-Ill trajectories 
contributing to certain final momenta around the caustic struc- 
ture are indicated. Particularly many type-Ill trajectories con- 
tribute at p c , forming the semi-classical analogue of the LES: 
a caustic. Analytically, the mapping {pox,Poz} -t {p x ,Pz} 
has a singularity at the caustic. In the numerical TC-SFA 
calculations this singularity manifests itself as an increase of 
the number of trajectories ending up close to p c . The same 
phenomenon is known for rescattering where the trajectories 
merge and branch at classical cut-offs [20]. Corresponding 
exact quantum mechanical calculations show an enhancement 
of the probability there fl9|,|2l|]. As is visible in Fig. [2^, panel 
III, the LES-caustic is a classical cut-off for type-Ill trajecto- 
ries with not too small asymptotic lateral momenta. Feature 
(ii) from above, i.e., the radial structures around p = 0, are 
due to interference of all four types of trajectories. The side- 
lobes forp z < (iii) emerge due to interference of type-II and 
type-Ill trajectories 12211 . Both features are not reproducible 
by purely classical ensemble simulations but require quantum 
approaches including Coulomb effects. 

In summary, we presented a semi-classical approach based 
on quantum orbits to calculate the strong field ionization ma- 
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